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ABSTRACT 



Aims. We investigate the linear stability properties of the plane interface separating two relativistic magnetized flows in relative 
motion. The two flows are governed by the (special) relativistic equations for a magnetized perfect gas in the infinite conductivity 
approximation. 

Methods. By adopting the vortex-sheet approximation, the relativistic magnetohydrodynamics equations are linearized around the 
equilibrium state and the corresponding dispersion relation is derived and discussed. The behavior of the configuration and the regimes 
of instability are investigated following the effects of four physical parameters, namely: the flow velocity, the relativistic and Alfvenic 
Mach numbers and the inclination of the wave vector on the plane of the interface. 

Results. From the numerical solution of the dispersion relation, we find in general two separate regions of instability, associated 
respectively with the slow and fast magnetosonic modes. Modes parallel to the flow velocity are destabilized only for sufficiently 
low magnetization. For the latter case, stabilization is attained, additionally, at sufficiently large relativistic velocities between the two 
flows in relative motion. 

Conclusions. The relevance of these results to the study of the stability of astrophysical jets is briefly commented. 
Key words. Kelvin-Helmholtz instability - relativistic MHD - plasma physics 



1. Introduction 

The Kelvin-Helmholtz instability (KHI henceforth), acting at the 
contact interface between two flows in relative motion, plays a 
dynamically important role in a number of different astrophys- 
ical scenarios such as, for example accretion disks, planetary 
magnetospheres, stellar and extragalactic jets. It is a classical 
instability in fluid dynamics, already discussed at the end of 
ninete enth century in th e work s of Von Helmholtz & Monats 
d 1868b and Lord Kelvin dl871l) . The classical results of the lin- 
ear analysis for incompressible flows separated by a planar vor- 
tex sheet, both in the presence of magnetic fields and for dif- 
ferent geometries are summ arized in Chandrasekhar's mono- 
graph (Chandrasek har 19611) . They have been later generalized 
for c ompressible case both in the pure hydrodynamical limit (see 
e.g. iGerwin 19681 : ISen 1964)) and in the magnetohydr odynamic 
(MHD) limit (see e.g. ISen 1961: iPu & Kivelson 1983b extend- 
ing the study to the supersonic regime. In particular it was shown 
that in the hydrodynamic case no unstable longitudinal (with re- 
spect to the plane interface) modes are allowed for flow Mach 
numbers larger than V8, while transverse modes are always un- 
stable. 

More recently Choudhury & Lovelace (1986) studied the lin- 
ear instability behavior in the case of a finite-thickness bound- 
ary layer, as a function of the Alfven velocity and inclination 
of the wave vector with respect to the interface. Their main re- 
sult was the proof of the existence of two distinct regimes of 
instability in the wavevector-Mach number plane correspond- 
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ing to standing and traveling MHD waves, the extent of the in- 
stability regions depending upon the magnetization parameter. 
Magnetohydrodynamic flows in c ylindr ic al and slab geometries 
have been studied by Ferrari et al. dl98ll) . dl982l) showing the ex- 
istence of a new range of unstable modes, arising for high Mach 
number flows, due to standing waves created by reflections at the 
boundaries of the structures. 

For astrophysical problems involving high-energy phenom- 
ena relativistic effects become important either because the rel- 
ative velocities are close to the speed of light or because the 
plasma is relativistically hot or because the magnetic field is very 
strong and the Alfven speed becomes close to the speed of light. 
For the non magnetic case a relativistic extension of the KHI 
study was considered, f or a pl ane parallel velocity discontinuity, 
by Turland & Scheuer (1 19761) . Blandford & Pringle d 19761) and 
more recently by Bodo et al. d2004l) . In particular Bodo et al. 
d2004l) showed that the dispersion relation can be solved analyt- 
ically in a frame of reference where the media in the two half- 
spaces have equal and opposite velocities and that the stability 
criteria have the same form of classical ones, simply substitut- 
ing the Mach number with the rel ativistic Mach number, that was 
first introduced by Konigl d 19801) . 

The relativistic MHD case still lacks of a full stability analy- 
sis and this is what we present in this paper. Early attempts were 
carried out by Ferrari et al. (1980). Their analysis, however, as- 
sumed the limit of Alfven velocity much smaller than the speed 
of light, thus allowing to neglect displacement currents, while 
the speed of the flow could reach relativistic values. The main 
result of this study was the confirmation that KHI disappears 
when the Alfven velocity approaches or becomes larger that the 
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speed of sound and when a large density contrast exists between 
the two media. 

Relativistic magnetized flows can be found in astrophysical 
environments associated to AGN jets, to jets from galactic X-ray 
binaries, such as GRS 1915+105 (Mirabel & Rodriguez 1999), 
and to pulsar wind nebulae (PWNe). Applications of the Kelvin- 
Helmholtz instability mechanism to the case of AGN and galac- 
tic microquasar jets is widely discussed in the literature, mainly 
for cy lindrical geometries, such as the recent analysis by Hardee 
d2007l) . who studied axially magnetized cylindrical relativistic 
jet embedded in a magnetized sheath. In addition KHI appear to 
be the physical mechanism to slow down relativistic jets b y pro - 
ducing entr ainmen t of the ambient medium. Bodo et al. ([2003), 
Rossi et al. (12004 and Rossi et al. (2008) have studied the non- 
linear evolution of KHI in cylindrical relativistic jets to explain 
the dichotomy of FRI - FRII extragalactic radiosources. Instead 
the planar case can apply to the case of PWNe, when a rela- 
tivistically hot magnetized plasma forms where the pulsar ultra- 
relativistic wind interacts with the surrounding supernova ejecta. 
The linear stability analysis that we present here will be relevant 
for understanding the physical phenomena involved in the desta- 
bilization process. 

In this work we present, for the first time to our knowledge, 
a complete treatment of the linear stability problem in the plane- 
parallel geometry, by considering the full system of equations 
(i.e. without neglecting the displacement current and account- 
ing for non relativistic as well as relativistic Alfven velocities). 
It is shown that the stability properties depend on four physical 
parameters which we conveniently choose to be the relativistic 
Mach number, the flow velocity, the Alfven speed and the incli- 
nation of the projection of the wave vector onto the plane of the 
interface. 

The paper is organized as follows. In {32.2l we derive the gen- 
eral dispersion relation, in |3]we present and discuss the results 
for the non-relativistic and relativistic cases, and we add some 
astrophysical comments in ^4] 



metric tensor. For ease of notation, factors like An have been set 
to unity in the derivation. 

The electromagnetic field tensor F 1 ™ appearing in the defini- 
tions of the stress energy tensor, Eq. ©, obeys Maxwell's equa- 
tions 



d[ a F, 



jJV] 



0, 



-r. 



(3) 



where [...] denotes anti-symmetrization and J v is the charge four- 
current. In the limit of infinite conductivity, the rest frame elec- 
tric field vanishes identically and the electromagnetic field tensor 
becomes orthogonal to the fluid four-velocity, i.e. F^up = 0. 
Note that this is the only approximation introduced in ideal 
RMHD and it becomes identical to the well known non- 
relativistic expression E = -vxB, where E is the electric vector 
field. Unlike classical MHD, however, the displacement current 
is not discarded in RMHD and explicitly enters through the def- 
inition of the current, given by the second of Eqs. (O. The high 
conductivity limit allows to write F 1 ™ = ^^b^Uy, where V 1 is 
the magnetic induction four- vector: 



\f = 



B 

yv ■ B, h y (y ■ B) v 

7 



(4) 



with B and e"^ denoting the rest frame magnetic field and the 
Levi-Civita pseudo-tensor, respectively. 

For the purpose of our derivation, we explicitly rewrite Eq. 
([TJ in components: 



dU ydf_ 



0. 



(5) 



where U and f are the vector of conserved variables and corre- 
sponding fluxes: 

U = [py, Wt y 2 v j - b°b j , w,y 2 - b°b° - p,] T , (6) 
/'' = [pyv\ w f yVV - bV + p t 6 iJ , w t y 2 v l - b% l f , (7) 



2. Relevant Equations 

2.1. The Equations of Relativistic MHD 

The equations of relativistic magnetohydrodynamics (RMHD 
henceforth) govern the evolution of a (special) relativistic mag- 
netized fluid. They can be cast as a system of conservation laws, 
describing energy-momentum and mass conservation: 



V^o, v^(p^) = o, 



(1) 



where p is rest mass density and u M = (y, yv) is the fluid four- 
velocity (y = Lorentz factor, v = three velocity) such that 
u^u^ - -1. The speed of light is set to unity everywhere. The 
expression of the energy-momentum tensor can be found under 
the physical assumptions of constant magnetic permeability and 
infinite conductivity, appropriate for a perfectly conducting fluid 
(see, f or example, the books by Anile ( 119891) and Lichnerowitz 
(fl967l» . 

Under these conditions, the stress energy tensor for a perfect 
fluid interacting with an electromagnetic field decomposes into 
rpfiv _ Tp^+TgTj, where the fluid (FL) and electromagnetic (EM) 
contributions are given, respectively, by 

T^=phu"u v +pTr, T^F^-^rfF^. (2) 

Here h and p are used to denote the gas specific enthalpy and the 
thermal pressure, while rf v = diag (-1, 1, 1, 1) is the Minkowski 



where p, and w, are, respectively the total pressure and enthalpy, 
expressed as the sum of thermal and magnetic contributions: 



Pi =P + ^\b\ 2 , w,=ph + \b\ 2 : 



(8) 



where \b\ 2 = ¥'b a . Magnetic field evolution is given by the in- 
duction equation in Maxwell's law: 



dB 

dt 



- V x (v x B) = , 



(9) 



which conserves the same form as in the classical case. Proper 
closure is provided by specifying an equation of state which we 
take as the constant F-law: 



h = 1 



r p 



(10) 



r-ip' 

where T is the polytropic index of the gas. 
2.2. Dispersion Relation 

Our setup consists of a planar vortex sheet interface in the xz 
plane at y = separating two uniform flows moving in the x 
direction with opposite velocities v(y > 0) = fti and v(y < 0) = 
-pi. The fluids have equal density and pressure and are threaded 
by a uniform longitudinal magnetic field along the direction of 
relative motion, B = Boi. 
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k± — y(k + to/3) , l± — l± , m± = 



We start our analysis by introducing small deviations around 
the equilibrium state. This is better achieved by working in the 
rest-frames of the fluids, where perturbations are sought in the 
form: 



3 

(25) 



so that a generic perturbation in our original frame can be ex- 
pressed as 



B = B +B + .. 
v = v' + ... , 
p =po +p' + ■■■■ 
p = Po + P + - 



(11) 

(12) 
(13) 
(14) 



oc exp [i (kx + l±y + mz — u>t)] . 



(26) 



where the tilde denotes quantities in the rest frame. The per- 
turbed terms of the first order can be expressed by the following: 



"T^ oc exp (k ± x + l ± y + m ± z - <i ± f)J 



(15) 



where W' ± = (B ± , v' ± ,p' ± , p' ± ). Subscripts + and - correspond to 
the regions y > and y < 0, respectively. Wave numbers and 
frequency are denoted with k±, l ± , m ± and u> ± . 

Proper linearization of the equations leads to the dispersion 
relation (Komissarov d!999l) ): 



u 2 ± fi ± 



(k 2 ± + Ii + ml) 2 k 2 ± + li + 



~ 2 
1)11 



= 0, 



P± 



k 2 +P 

= c 2 + v 2 A -c 2 vl- 1: = 



%± + l± + m ± ' 



v± 



c 2 vU 



C s = 



Va = 




Jpoho + 



which allow the propagation of the Alfven mode 



(16) 
(17) 
(18) 

(19) 
(20) 



By substituting Eqs. ( l24b ) and (l25T l into Eq. (f2Tb we derive 
immediately that the possible roots are always real, hence we 
may already conclude that the Alfven mode will not contribute 
to the instability. In the case of Eqs. d22l i and (l23l we see instead 
that the resulting slow and fast magnetosonic modes include both 
real and imaginary roots, and therefore we expect instability in 
the parameter domains corresponding to these modes. 

Direct substitution into the induction equation (0 and into 
the y and z components of the momentum equation in (|5]l yields 

Bo— B' y± + [a)B 2 + ypohou*] v' y± - l ± (p' ± + B B' X± ) = , (27) 



Boy B' z± + [a>B 2 + yp h a>*\ v' z± - m(p' ± + B B' X± ) = , (28) 

(cj + kB)B' x± - B (l ± v' y± + mv' z± ) = , (29) 

(w + kj3)B' y± + kB Q v' y± = , (30) 

(oj + kp)B' z± + kB Q v' z± = , (31) 

where k T and co+ are given by Eqs. (I24ll25l l. 

With the aid of Eq. (0 one may express the total pressure 
perturbation as: 



P±=P± + 



l±B 2 v[. ± 
to + kB' 



(32) 



and imposing the displacement matching condition at the inter- 
face: 



y+ 



kB 



kB 



and the total pressure continuity {p' + 
obtains (together with Eq. |27) 

l ± _ pp/ipy V - kB) 2 + (gj 2 - k 2 )B 2 
and the two magneto-acoustic modes, (i) slow magnetosonic L pohoy 2 (tL> + kB) 2 + (oj 2 - k 2 )B^ 



k 2 A 

K A± 



(21) 



(33) 

p'_) at the interface, one 
(34) 



kl + n ± + m 2 s± 2— 2 
and (ii) fast magnetosonic 



1/2 



(b 2 f± 



kl + ^+ml 



1/2 



(22) 



(23) 



Introducing the dispersion relation of slow and fast magne- 
tosonic waves (see the right bracket of Eq. ( fT6l l) and solving for 
I 2 one gets 



-m 2 + 



co 2 



l[(C 2 + V 2 A )kl-col]-C 2 V 2 A kt 



C 2 V 2 .(kl + u>l) - (C 2 + V 2 A )cbl 



(35) 



Here C s and Va are the sound speed and the Alfven speed respec- 
tively. Notice that both C s and V A are invariant under Lorentz 
boosts in the x-direction. 

In order to obtain the dispersion relation in the laboratory 
frame, one has to transform all the quantities by mean of a 
Lorentz transformation: 



oj ± = y(u> + kB) , 



(24) 



which, combined with Eq. d34b . gives the desired dispersion re- 
lation of modes connecting through the planar interface. 
A number of considerations can be immediately drawn. 

- The dispersion relation obtained by substitution of Eq. ( f34b 
into Eq. ( f35l > is an 8-th degree polynomial in u>. However, 
only 4 out of a total of 8 complex roots do actually satisfy 
Eq. d34l l and thus are physically significant. 
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Roots with positive imaginary part of a> identify unstable 
modes. The growth rate will be more conveniently expressed 
by the imaginary part of the dimensionless quantity 



cf> = 



(36) 



Bq and 



In the non-relativistic limit (i.e. p — > 0, ph » p 
h — > 1) our dispersion rel ation r educes to the form given 
by Choudhury & Lovelace ( 1986) in their equations (7) and 
(A2a); 

In the limit of vanishing magnetic field Eqs. d34l and d35l > 
simplify considerably and the dispersion relation reduces to 
the one derived by Bodo et al. (2004) (see their Eqs. (l33l > and 
03); 

For the general case the system of equations Eqs. ( 1341 [35l > 
has to be solved numerically and these results are discussed 
in the following sections. 



3. Results 

Our study will examine the instability dependence on 4 param- 
eters: 1) the fluid velocity (3 (alternatively, the Lorentz factor y), 
2) the relativistic Mach number 



Mr = 



(37) 



3) the Alfvenic Mach number £ = V A /C S and 4) the ratio / = 
m/k, that gives the angle between the wave number projection on 
the xz plane and the flow velocity. In general, relativistic effects 
come into play whenever one of /3, M r or ( (or a combination 
of them) describes situations of high fluid velocities, hot gas or 
strong magnetic field, respectively. To this purpose, it is useful 
to use the explicit relation linking the Alfven speed and the other 
parameters: 



V A = 



ft 



Vy8 2 +M r 2 d-y6 2 ) 



(38) 



In the following, we will discuss the non-relativistic case and 
subsequently the relativistic case, both for oblique propagation 
(/ + 0) and for propagation parallel to the velocity interface 
(/ = 0). The latter case will be useful for the comparison with 
the results obtained by previous investigators. 

3.1. Non-relativistic flows 

As a starting point, let us now consider the case of a non- 
relativistic fluid, with /3 = 10~ 3 and r = 5/3. 

A convenient way to gain a comprehensive view of the 
general behavior of the instability, when one a llows for non - 
parallel propagation (i.e. f + , see also Gerwin dGerwin 1 968), 
Choudhury & Lovelace (1986)), is to represent the results on a 
contour plot such as the one in Fig.Q] where the contour levels 
show the growth rate Im(</>) as a function of / and M t , for fixed 
Alfven Mach number if = 1.6 (top panel). From the plot, one 
derives that i) two separate instability regions are present in the 
plane (A4 r , /), corresponding to slow and fast magnetosonic un- 
stable modes, and ii) they do not extend to the horizontal axis, 
i.e. at / = 0. By expanding the solution of the dispersion relation 
in terms of 1/f, one can show from Eqs. ([34] [35} that the cor- 
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2.0 2.2 2.4 
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Fig. 1. Two contour plots are shown. Top panel: the contour plot 
of constant lm(0) surfaces as a function of the inclination fac- 
tor / and relativistic Mach number Ai r for fixed Alfvenic Mach 
number if = 1.6 and flow velocity /3 = 10~ 3 . Bottom panel: the 
contour plot of constant Im(</>) surfaces as a function of the in- 
clination factor / and j3 for the fixed Alfven speed (Va = 10" 3 ) 
and the sound speed (C s = 6.3 x 10 4 ). Here F = 5/3 is used. 

responding first order term of the growth rate can be expressed 
by: 



v 2 A -p 2 



1 



(39) 



Eq. (1391 demonstrates that, for large values of /, the instabil- 
ity appears when Va < P- In Fig. Q~|(bottom panel) we show the 
contour plot of constant Im(</>) surfaces as a function of / and /? 
for Va = 10 3 and C s = 6.3 x 10 4 (this value corresponds to 
Air = 1 .6 (top panel), see Eq. [37l i. and the figure shows that the 
instability appears for /3 > 10~ 3 , i.e. when ft exceeds the value 
of Alfven speed parameter. We note as well that, while the re- 
sult in Eq. (l39l has a general validity, in the non relativistic limit 
/3 <s 1, the condition for instability derived above for large f's 
reads also f < Ai r , as evident for Fig.Q](top panel). 

For smaller f's, i.e. for small inclinations of the propagation 
direction with respect to the velocity interface, Fig.[T]shows that 
for / < 0.2 the system becomes stable for the set of parameters 
considered. 

This general behavior is better understood by looking at the 
real and imaginary parts of the roots giving rise to unstable 
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1,2 1.4 1.6 1.8 2.0 2.2 2.4 



Fig. 2. Imaginary (dotted-dashed lines) and real (solid lines) 
parts of the roots satisfying the dispersion relation as functions of 
the relativistic Mach number. Only roots associated with unsta- 
ble modes are plotted. For the sake of clarity, unphysical roots 
have been plotted as well. The other parameter set: (3 = 10~ 3 , 
£ = 1.6, / = 1 and T = 5/3. Instability arises when a < M r < b 
and Mr > c. 
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Fig. 3. Same as Fig. [1] but for ( = 1 . 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Fig. 4. Dependence of lm(<p) on the relativistic Mach number 
for fixed (3 = 1(T 3 and / = (T = 5/3). The different labels 
above each curve give the values of the Alfvenic Mach number 
<T= [0,0.2,0.5,0.8,0.9]. 
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Fig. 5. Instability Mach range AvVt r as a function of the magnetic 
field in the case of parallel propagation for ft = 10~ 3 (solid line, 
r = 5/3) and f3 = 0.7 (dashed line, F = 4/3). For the sake 
of comparison, the latter has been reduced by a factor of 50. 
Clearly, the instability is suppressed when the Alfvenic Mach 
number equals to 1 (for^ = 10~ 3 ) and 0.35 (for/? = 0.7). 



modes, shown in Fig. [2] This plot corresponds to a horizontal 
cut across the contour lines of Fig. Q]for fixed / = 1. Indeed, 
for M t < a, all roots are real yielding a stable configuration. For 
a < M r < b, however, the roots become complex conjugate with 
nonzero imaginary parts, making this range unstable. Between 
b < M r < c we have again a stable configuration, whereas for 
Mr > c the roots become purely imaginary and a second unsta- 
ble region can be recognized. It can be verified by direct substi- 
tution into Eq. ( l22b and Eq. ( 1231 ) that the two family of solutions 
separately correspond to the onset of slow (for a < M r < b) and 
fast {Mr > c) magnetosonic modes, respectively. 

For decreasing magnetization (i.e. lower values of (), the 
two instability regions approach each other, gradually shifting 
towards lower values of Mr- At the same time, lower values of / 
become prone to instability, see Fig. [3] The limiting case / = 
is reached for f < 1, when the previously identified instability 
ranges merge into a single unstable region. 



We will now examine the case of parallel propagation (/ = 
0) for different values of £ < 1. These results may be directly 
compared with those obtained by Ferrari et al. ([1980). In Fig. |4] 
we plot the growth rate as a function of Mr for several values of 
the Alfven Mach number in the range < (, < 0.9. For vanishing 
magnetic field (£ = 0) the instability disappears when M r > V2, 
consistently with the results shown by Bodo et al. ((2004))). As 
the magnetic field increases, higher values of £ induce a stabi- 
lizing effect by raising the pressure and forcing the flow to be 
channeled along the field lines. This has two noticeable effects: 
one is to decrease the maximum value attained by the growth 
rate and the other is narrowing down the instability Mach range 
AM, = M™ x - Mf n , where A4™ x (A1, min ) is the largest (small- 
est) value of the Mach number above (below) which the config- 
uration is stable. Note that the values of the sound speed at the 
growth rate maxima are not relativistic. 
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Fig. 6. Top panel: contour plots of constant lm((p) lines in the 
Air - f plane for (3 = 0.7, ( = 1.6. The two instabilities branches 
are clearly visible. Bottom panel: contour plots of constant lm(0) 
lines in the - f plane for V A = 0.7, C s = 0.43. Here T = 4/3 
was used. 



Fig. [5] shows the dependence of the critical range of the rel- 
ativistic Mach number on the magnetic field (solid line). As it is 
clear from the plot, when ( — > 1 , the relativistic Mach number 
range tends to zero and the instability disappears. For f > 1 the 
interface is always stable. Our res ults are thus in full agreement 
with those of Ferrari et al. (11980b for the same parameter range 
(J3 = 10- 3 , / = 0). 

3.2. Relativistic flows 

In the previous subsection all the relevant quantities were de- 
scribing essentially flows in non-relativistic conditions. We now 
consider situations in which the fluid exhibits relativistic behav- 
ior, at least in some parameter range. 

The top panel in Fig.[6]shows the contour levels of the growth 
rate for fixed flow velocity /3 = 0.7 and Alfvenic Mach number 
£ = 1.6. We have again a slow magnetosonic modes instability 
region for 1.63 < vVt r < 1.8, and a fast magnetosonic modes 
instability one for Ai r > 2.02. For large /'s the latter one only 
survives and the discussion for the non relativistic case is still 
valid, i.e. the instability condition is Va < P (see Fig. [6] bot- 
tom panel). In fact, a direct comparison between Figs. |6]and[T] 
reveals a similar qualitative behavior between the two instabil- 



Fig. 7. Contour plots of constant Im(</>) lines in the Ai r - f plane 
for ft = 0.7, £ = 1.2 (top panel) and f = 1 (bottom panel). Here 
F = 4/3 was used. 



ity regions corresponding to slow and fast magnetosonic modes. 
Indeed, in the relativistic case, higher values of Mi are necessary 
to promote instability and the gap separating the two regions be- 
comes larger. However, for lower values of only one instabil- 
ity region survives and there are no contour lines for Ai r < 1 .4, 
see the top panel in Fig. [7] This is a consequence of the fact 
that higher values of the relativistic Mach number correspond to 
lower sound speeds and, since the latter cannot exceed the up- 
per limiting value of 1 / V3, M r has a lower physical cutoff at 
V2/?y 1.386. Thus, no merging can take place as f decreases, 
since the leftmost region disappears below this threshold. This 
effect is clearly visible in Fig. [7j In the limit of parallel propa- 
gation (/ = 0), the growth rate decreases with stronger magne- 
tization, see Fig. [8] This is also manifested by a reduced range 
of the Mach number values for which instability exists. This be- 
havior has already been discussed for non-relativistic flows in 
{sl3.ll where the flow was shown to become stable for values of 
the Alfven speed close to the speed of sound. For the present 
case (yS = 0.7), however, stability is approached when the Alfven 
speed Va — 0.35C S , see Fig. [8] This tendency can also be recog- 
nized by the profile of AAi r (the critical Mach number instability 
range) as function of £ shown in Fig. [5] (dashed line). As antic- 
ipated, the instability is quenched for ( =t 0.35, an effect which 
can also be justified by the increased kinematic inertia. 
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Fig. 8. Effects of magnetic field on the growth rate. Labels on 
each curve indicate values of (. The set of parameters is: p — 0.7, 
f = [0, 0.2, 0.3, 0.34], / = and T = 4/3. For f = 0.35 the flow 
is completely stabilized. 




0.62 0.64 0.66 0.68 0.70 0.72 



Fig. 9. Effect of the kinematic factor on the instability when 
/ = (T = 4/3). As it is clear, for p 0.71 the instability is 
suppressed by relativistic effects independently on M T . 



Th ese results confirm the general trend already discussed in 
d 1 980): however differences in the growth rates are found when 
the Alfven velocity approaches the velocity of light and displace- 
ments currents become relevant. 

In Fig. [9] indeed, we show the critical value of f (above 
which the instability is suppressed) as a function of the flow ve- 
locity p. Thus £ cr monotonically decreases from ~ 1 (at/J = 0.6) 
to for p =i 0.71. Higher value of p leads to a stable interface, 
even without magnetic fields. This is a major difference from the 
classical MH D case . The latter result was already introduced by 
Bodo et al. ((2004)), where it was shown that positive growth 
rates are subject to the condition 



P<Q 



(40) 



Since the right hand side is a monotonically increasing function 
of C s , we conclude that for p > 0.7071 (C s = 1/ V3) the insta- 
bility is suppressed by kinematic effects only, regardless of the 
value of the magnetic field. 



4. Summary 

We carried out a linear stability analysis of the relativistic mag- 
netized Kelvin-Helmholtz instability problem in the vortex sheet 
approximation. Solutions to the dispersion relation have been 
sought in terms of four parameters giving the strength of the 
magnetic field, the relativistic Mach number, the flow velocity 
and the spatial orientation of the wave vector. The main results 
of this analysis can be summarized as follows: 

1 . We have examined the cases of non-relativistic and relativis- 
tic flows separately. We have found that, for fixed values 
of the Alfven velocity, two separate regions of instability 
appear in the plane (M T ,f). These two instability regions 
are associated with the destabilization of slow and fast mag- 
netosonic modes. For large /'s, only the fast magnetosonic 
modes are unstable in the region Va < P- 

2. For high values of the Alfven velocity, modes propagating 
parallel to the flow velocity are stable. As this velocity de- 
creases, the two instability regions corresponding to slow 
and fast magnetosonic modes tend to merge and the parallel 
modes eventually destabilize. This general behavior is com- 
mon for both non-relativistic and relativistic flows. 

3. We have found that slow modes non-parallel to the flow di- 
rection are gradually excluded from the instability plane as 
they become unphysical. This effect takes place for increas- 
ing relativistic velocities or low magnetic fields or a combi- 
nation of both. 

4. When considering the case of parallel propagation (/ = 0), 
we have found that, similar to the non-magnetic counterpart, 
the flow becomes linearly stable when the relativistic Mach 
number exceeds a critical value. In the limit of vanishing 
Lorentz force, this threshold reaches the maximum value of 
V2, in a frame where the fluids have equal and opposite ve- 
locities. 

5. For increasing magnetic field strength, the maximum unsta- 
ble growth rate decreases and the instability exists for a nar- 
row range of Mach number values. 

6. The main difference from the non relativistic MHD is that, 
at higher flow velocities (for / = 0), kinematic effects stabi- 
lize the flow even for smaller values of the Alfven velocity. 
Furthermore, the computed growth rates attain lower values 
than their classical counterparts. When the flow velocity be- 
comes higher than 1 / y2, the flow does not need any more 
the magnetic field for stabilization. 

Relativistic flows stabilize the KH modes due to the concurrence 
of two effects: first, at high relative velocities, the mode cou- 
pling between the two half-spaces loosens because of causality 
effects, this effect is present in the non-relativistic case as well 
but is even more relevant at relativistic velocities; second, in the 
relativistic regime the inertia of the fluid particles is dynamically 
augmented, hampering the response of the fluid to the instability. 

Destabilization of relativistic Kelvin-Helmholtz modes has 
been applied in the astrophysical context of jets of extended ex- 
tragalactic radiosources and relativistic galactic sources to ex- 
plain the generation of observed morphologies and the driving 
mechanisms for their radiation emission. While originally these 
modes were thought to yield disruption of the collimated propa- 
gation, it was later realized that the head of jets creates a bow 
shock and a cocoon that shield the flow from steep gradient 
boundary layers while leading to generation of a turbulent state. 
Turbulence is functional in producing mixing of jet and ambi- 
ent matter with entrainment, angular momentum transport and 
slowing down of the flow (Rossi et al. 2008). The behavior of 
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KHI as a function of the physical parameters is likely to the at 
the base of the different classes of radiosources. In addition long 
wavelength modes may develop and generate global structures 
as the wiggles and knots that are observed in astrophysical jets. 
The present fully relativistic analysis confirms that in order to 
activate KH modes small Va/Cs ratios are needed. This condi- 
tion can be realized in particular by very hot plasmas, even for 
relatively large magnetic fields, as it appears to be the case in 
high-energy sources. 

KHI appears also to be relevant in the case of PWNe; in fact 
in these objects a relativistically hot magnetized plasma is pro- 
duced by the pulsar ultra-relativistic wind and interacts with the 
surrounding supernova ejecta. These PWNe are internally struc- 
tured and presents high velocity flow channels that can become 
KH un stable. As already mentioned, Bucciantini & Del Zanna 
(2006) have examined numerically the nonlinear development 
of KHI in these channels, assuming planar geometry, for inter- 
preting the signature of the instability on the synchrotron emis- 
sion. The present linear stability analysis is complementary to 
their study as it allows a better comprehension of the physical 
phenomena involved in the destabilization process. Again the 
presence of a hot plasma appears necessary to start the channel 
destabilization. 
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